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Abstract 


In the recent hydroelectric power plants built in Mexico (P.H. El Cajón, 
P.H. La Yesca, and Grijalva River tunnels), the diversion tunnels have 
been designed with a combination of materials: hydraulic concrete in its 
template and concrete cast in ¡ts walls and vault, which allows to reduce 
hydraulic losses and get a greater discharge capacity to guarantee the 
hydrological safety of the dam during the build time, in addition to 
achieving considerable cost savings. To monitor this capacity, several 
instrumentation projects have been carried out, recording the hydraulic 
variablesin different sections, in addition to estimating the representative 
parameter of the composite roughness. 


This article relates how estimates the composite roughness 
parameter in diversion tunnels, in a case study: "Communication tunnels 
of the Grijalva River”, upstream of the Penitas Hydroelectric Power Plant, 
using hydraulic instrumentation to measure, among other variables, the 
speed in the area near to the wall. And applying the boundary layer 
theories of Nikuradse and Prandtl-Von Kármán to determine the 
equivalent roughness in the walls, as well as the review of 17 empirical 
equations to finally obtain the representative composite roughness of the 


tunnels. 
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Resumen 


En los recientes proyectos hidroeléctricos construidos en México (P.H. El 
Cajón, P.H. La Yesca y túneles del río Grijalva) se han diseñado las obras 
de desvío con una combinación de materiales: concreto hidráulico en su 
plantilla y concreto lanzado en sus paredes y bóveda, lo que permite 
reducir las pérdidas hidráulicas y obtener una mayor capacidad de 
descarga para garantizar la seguridad hidrológica de la presa durante el 
tiempo de construcción, además de lograr considerables ahorros en los 
costos. Para monitorear dicha capacidad se han realizado diversos 
proyectos de instrumentación, lo que permite registrar las variables 
hidráulicas en diferentes secciones, además de estimar el parámetro 


representativo de la rugosidad compuesta. 


Este artículo trata sobre la estimación del parámetro de rugosidad 
compuesta en túneles de desvío con caso de estudio: “túneles de 
comunicación del río Grijalva”, localizados aguas arriba de la Central 
Hidroeléctrica Peñitas, utilizando instrumentación hidráulica que permite 
medir, entre otras variables, la velocidad en la zona cercana a la pared. 
Y aplicando las teorías de capa límite de Nikuradse y de Prandtl-Von 
Kármán determinar la rugosidad equivalente en las paredes, así como la 
revisión de 17 ecuaciones empíricas para obtener finalmente la rugosidad 
compuesta representativa de los túneles. 
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Introduction 


The Federal Electricity Commission (CFE) through the Hydropower 
Projects Coordination (CPH) along with the Hydraulics Laboratory have 
carried out several hydraulic instrumentation projects since 2008: in the 
Cajón dam, Yesca dam, and Grijalva River tunnels. In the Yesca dam 
several structures were instrumented to study their hydraulic behavior: 
the diversion tunnels; the generation elements, the discharge channel and 
oscillation chamber; the spillway; as well as the river downstream of the 
dam; and the reservoir. These instrumentation projects and studies 
consist mainly of the installation of pressure sensors or cells and Doppler 
gauges to calculate velocity and flow rate, as well as other variables such 
as the study of the composite roughness. 
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In 2007, as a result of the 1077 mm of rainfall accumulated from 
October 28th to November 3rd in the Grijalva River in Chiapas, there was 
a landslide of 55 million cubic meters of earth and rock on the right bank 
that caused a blockage in the river 800 m long and 300 m wide, posing a 
serious safety risk forthe Penitas dam which in the case of breakage could 
generate an unregulated flood. In addition to this, serious damage to 
downstream and upstream populations could be caused due to the 
reservoir overflow. Moreover, afterthe landslide, the connection between 
the reservoirs upstream of the fall and the Penitas dam was reduced to a 
channel measuring 930 m long and 70 m wide (Figure 1). 


Estación H. Juan de Grijalva 
sj (Vaso Superior) 


Portal de 
Entrada 

Ko 0+105 
*.(RISONIC) 


Túnel 


Canal 1 


Juan de Grijalva 


Portal de 
Salida 


Estación Hidrométrica 
(Vaso Inferior) 


Figure 1. Location of the communication tunnels of the Grijalva River 
and the Juan de Grijalva channel afterthe landslide. 
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Among the actions taken to deal with flows of up to 3 500 mY/s, 
were: 1. Building a channel in the reduction zone along the course of the 
river; 2. A temporary suspension of the generation of the Malpaso 
(upstream); 3. Building two 14x14 m portal section tunnels to 
complement the capacity of the river between the Malpaso and Penitas 
dams. The 3rd option for the tunnels requires monitoring for efficiency in 
hydraulic conduction, as well as in natural geological movements of 


structure, through hydraulic and geotechnical measurement instruments. 


Hydraulic instrumentation was used to measure hydraulic levels, 
flow rates, and velocity profiles in defined sections, both in low water and 
rainy periods. The composite roughness parameter representative of the 
tunnel was also studied, which would allow for the hydraulic capacity to 
be defined. 


The main objective of this work is to present a methodology to 
obtain the composite roughness parameter representative of the tunnels 
in the Grijalva River from direct measurements near the wall and apply 
the Prandtl-Von Kármán boundary layer theories and the Nikuradse 


equivalent roughness. 


The specific objectives are 1) to study the hydraulic performance of 
the tunnels using hydraulic instrumentation; 2) to estimate the wall 
roughness parameter from velocity measurements in the boundary layer 
zone; 3) to estimate the composite roughness parameter through point 
measurements and applying Prandtl-Von Kármán and  Nikuradse 


boundary layer theories; 4) to analyze a mathematical model and the 
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minimum square error (MSE) to find the best empirical equation for the 
composite roughness calculation, and 5) to obtain a roughness-flow 


curve. 


This will aid in the design of hydraulic tunnels with composite 
roughness as well as help it in its efficiency with discharge hydraulic 
capacity. This research was carried out under real operating conditions. 


Background 


Safety of tunnels and the composite roughness 


According to Marengo (2019), some considerations must be taken into 
account when designing and constructingtunnels. There is the experience 
during the construction of the Aguamilpa Hydropower Project: two 
extraordinary floods occurred and put the design conditions of the dam 
and the diversion tunnels to the test. It consisted of two 16 x 16 m portal 
section tunnels, a 1 100 m and 1 200 m long (Tunnel 1 and Tunnel 2, 
respectively), for a maximum flood of 6 700 m/s, in additionto a 10 m 


fusible tunnel and a channel on the right bank. In 1990, from August 11 
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to 22, there was a flood of 5 300 m3/s which caused an increase in the 
reservoir of 102.60 m above sea level. In 1992, from January 16 to 20, 
there was another flood of 10,800 m/s of instantaneous discharge and a 
maximum daily discharge of 9 334 m3/s. This caused the reservoir level 
to rise from 70 m above sea level to 108 m above sea level, which 
provoked the upstream overflow. The tunnel fusegate had to be opened 
and flooding began between the cofferdam and the dam to prevent further 
damage and the collapse of the structure. Therefore, the maximum level 
of 123.60 meters above sea level (MASL) was reached, just below the 
face of the dam. After that, there was a second flood of 7 700 m/s, 
reaching a height of 112.40 MASL, filling the enclosure between the dam, 
which caused a delay in the construction program of the Dam. 


The aforementioned events left the following lessons and 
experiences in regards to the tunnels: 1) the importance of placing a 
fusible tunnel to safeguard the integrity of the dam; 2) the advantage of 
using hydraulic concrete only in the tunnels' floor which allows the return 
period to change from 164 to 188 years (Tr). Likewise, by changing to 
concrete in the walls and vault, the Tr can increase to 313 years, which 
represents a significantincrease inthe safety of the hydraulic, geological, 
and structural system, also allowing for significant savings. Out of this, 
the concept of "composite roughness" is born that can be implemented 
with a risk analysis for each instance of the tunnels' design. 
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General aspects of composite roughness in tunnels or 


channels 


No matterthe shape of the subsections of the cross-section of a channel 
or tunnel (circular, trapezoidal, arch, horseshoe, etc.) when it comes to 
composite roughness, the main factor is the resistance coefficient n, 
which can be combined in the floor, walls, and vault. This coefficient 
modifies the velocity distribution along the cross-section (Marengo, 
2019), and for its study, several equations allow for the definition of its 
value homogeneously considering only the wetted perimeter (Marengo, € 
Arreguín, 2008). The study of the effects of varying the roughness on flow 
has been estimated through several types of analyses: physical, 
numerical, 2D, and three-dimensional models, as well as by on-site 
instrumentation combined with theoretical and empirical numerical 


analyses, which is the subject of study of the present research. 


Flow resistance depends on four components (Rouse, 1965): 1) the 
friction layer, 2) drag, 3) deformation resistance and 4) resistance 


associated with flow instability. Then Weisbach's drag coefficient f, is a 


function (F) of: 


f =F(Re,k,n, N, F, U) (1) 
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Where: Re is the Reynolds number; k is the relative roughness, 
which is expressed as ks/Rhk, where ks is the equivalent of the wall 
surface roughness and Rh is the hydraulic radius; y is a function of the 
geometric shape of the cross-section; N is the uniformity number of the 
channel for both the contour and the plane; F' is the Froude number, and 
U is the degree of instability of the flow. Rouse (1965) showed that the 
Moody diagram is a special case of Equation (1) for steady flow in rigid 
straight pipes of constant diameter by only considering two of the six 
independent parameters, which are: the Reynolds number and the 
relative roughness ks/Re. The latter allows the calculation of Nikuradse 


roughness. 


The formulas commonly used in open channels with a uniform flow 


to calculate the resistance coefficient are: 


Manning's equation: 
V == rpBs; (2) 
n 
Homogeneous Manning's equation: 


1/2 
Y = 22 pp2B si (3) 


Mg 
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Darcy-Weisbach Equation: 

ve aupa D1/2 (pm (4) 
ñ Ñ 
Chézy equation: 

P=CRAS (5) 


Hazen Williams equation: 


V = KuwChwR** S79* (6) 


Where V, is the average cross-sectional velocity; n, f, and C are the 
Manning, Weisbach, and Chézy drag coefficients, respectively; R, is the 
hydraulic radius; S¿ is the frictional slope; g is the acceleration of 


gravity. K, = 1m3/?/s. 


And the resistance coefficients can be related using the above 
equations: 


O A (7) 
f Vg Vg n ng [gRhSf 
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According to Marengo (2019), in professional practice, the Manning 
equation is still used. In the case of a uniform flow and from the spatial 
point of view (Chow, 1959), the flow, hydraulic area, velocity and 
discharge in each section of the tunnel are constant, and the hydrostatic 
slope S;, water surface S,,, and bottom of the channel S, are parallel, ¡.e.: 


Sw = Sp = Sy. And the denominatoris called the shear velocity, V,, which is 


defined as: 


This allows the coefficient of friction Cr to be determined: 
yy? 
C = () (2) 


According to Equation (9) and Equation (7), the following can be 


expressed: 
f_9_9g ve no 
f 8 c? K? dd el ( ) 
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When considering various roughness materials, for example, 
hydraulic concrete in floor and walls and shotcrete in the vault, or rock, 
the value of n, from Equation (2), an equivalent coefficient must be 


obtained and this could change with the water level (Marengo, 2019). 
Experimental design 


Instrumentation for measuring hydraulic variables 


Geometric and hydraulic characteristics of tunnels 


The research is centered on the case study of the communication tunnels 


of the upper and lower Grijalva River, upstream ofthe Penitas Dam (Angel 
Albino Corzo), with the characteristics mentioned below. 


It consists of two tunnels (Ti and T2) located on the left bank of 


the Grijalva River, with an arch section of 14 m wide by 14 m high and a 
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horizontal hydraulic concrete floor and walls, and a vault of shotcrete 
(Table 1). Both channels contribute to the regulation of the operation 


between the two dams Malpaso and Penitas (Table 2). 


Table 1. Geometric characteristics of the tunnels. 


Table 2. Operating Levels of Penitas and Malpaso Dams. 


NAME Elevation 95.50 m 


NAMO Elevation 87.40 m 
NAMINO Elevation 85.00 m 
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Superior elevation of spillway gates 91.14 m 
Flow rate by the generation unit 349.00 mY/s 


Gcdo 2IOOAAAA SONS 


Instrumentation installation 


A hydraulicinstrumentation project was carried out consisting of installing 
53 hydrostatic pressure sensors (measuring cells), programmed to 
measure water levels. Each tunnel has 24 sensors spaced out at a 
distance of 50 m of one another (in each tunnel), as well as 5 additional 


ones along each of the entrance and exit portals (Figure 2). 
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Figure 2. Location of cells in the Tunnels (Tunnel 1 above and Tunnel 2 


below in the Figure). 


The sensors are of the "“vibrating wire” type, which isinstalledinthe 
concave part of the concrete that joins the wall with the template in the 
cross-section of the tunnel, this area is known as "“chamfer” (Figure 3). 
The sensors contain a vibrating wire or cable that detects the hydraulic 
pressure using a diaphragm, which transforms the resonance frequency 
into a hydraulic pressure value and in turn, into a water level value. The 
sensors are previously calibrated in the laboratory and the programming 
is done in the SCADA system to transform it into values of hydrostatic 


pressure or hydraulic levels in meters above sea level (MASL). 
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Concreto lanzado — 15 
con fibra metálica A 
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para conducción de cables 
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<] — Celda de Presión 
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Figure 3. Placement of cells within the "chamfer”. 


For flow water, one module of Doppler equipment of the "Time in 
Traverse (T.T.)” type was installed per tunnel on the walls. Each module 
is located between sections 0+091 and 0+105 (Figure 4), which kept the 
local effect of the inlet from having an influence. The module consists of 
four pairs (eight sensors) measuring at different heights: 2.5 m, 5 m, 6.5 
m, and 8 m (Figure 5). It works for a velocity range of +-20 m/s and 
channels with a width of 0.75 m to 100 m and a minimum to a maximum 
range of temperatures varying from -20 *C to +70 *C and with an 
accuracy variation of only 0.5 %. 
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Figure 4. Location of “Time in traverse” meters (in the plant). 
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Figure 5. Location of “Time in Traverse” gauges (in profile or section). 


The T.T. concept is because the equipment sends a sound beam to 
another instrument installed at 45 degrees at the same elevation on the 
opposite wall, measuring in each pair the time it takes to return from 
downstream to upstream, thus calculating the fluid velocity. Water has 
small amounts of particles so the operation of T.T. meters is favored. 


To measure velocities, close to the wall, four doppler modules were 
installed. They have a higher frequency (3000 kHz), which allows them to 
measure velocities from 0.1 mto5.0 m distance from the wall. Each piece 


of equipment has two sensors at 25%, which gives redundancy to the 
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measurement. They were placed in three sections of each tunnel: 
KO+050, K0+655, and K1+115 (Figure 6), and at different heights: 2.5, 
5, 6.5, and 8 m (Figure 7), to measure a general spectrum in the section. 
Each piece of equipment was programmed to measure from 10 cm to 1.5 
m away from the wall because the objective was to measure in the zone 
of influence of the wall, this had an accuracy variation of only 0.1 % and 
measures at depths of up to 30 m. Its velocity ranges up to 6 m/s. 


K0+050,0 


El 
3 


5/3) K0+555,0 
395/5K1+115,0 


¡lll 


Túnel 1 


Y 

A 
a | 
IS 
5] 
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5 
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K1+115,035)5 


ko+050,0 535 


Figure 6. Location of Doppler velocity meters (in the plant). 
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Figure 7. Location of Doppler velocity meters (in section). 


Gauging campaign at the inlet portal 


To verify the certainty of the hydraulic variables, a gauging campaign was 
carried out in the inlet portal of the tunnels (Figure 1). It was compared 
with the flow rate of the T.T. meters installedinside thetunnel. A Doppler 


moving profiler (ADCP) was used for gauging, it had the following 
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characteristics: a depth varying from 75 cm to 30 m; an accuracy 
variation of 0.5 %; velocities of +-20m/s; and a temperature range of -5 
“Cto +45 *C. 


Flow water rate approximations were obtained from the ADCP 
movable Doppler for the T.T. meter ranging from -3 to 8 %. This 
difference is acceptable upon adding the tolerances of both pieces of 
equipment. Therefore, it was concluded that the T.T. meter installed 
inside the tunnels gives a better quality of instantaneous measurement in 
each reading: its configuration is designed to adjust with the other pairs 
of meters in case one of them is missing; it has a 0.5 % absolute error 
margin. Also, hydraulic parameters obey the fundamental theories and 
laws of fluid mechanics applied to the distribution of velocitiesin a rough 
bottom channel. 


Methodology 


Part 1. A theoretical approach to the calculation of the 


Manning's composite roughness, n« 
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Once the quality of the flow and the velocity measurement data has been 
verified, we proceed with the theoretical approach for the calculation of 
the composite roughness nc, starting with Bernoulli's Equation and 
Manning's Equation (Sotelo, 2002), where we know all the variables 


except for the roughness (Figure 8). 


Cálculo 


Y Perfil tipo H2 


Figure 8. Hydraulic calculation approach for roughness estimation. 
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If we only do not know the friction losses, hr: 
_ 2 pe 2 qn = AX 2nm2 1 1 
SL A 


Where S;, is the friction slope; 4, is the hydraulic area; R,,, is the 
hydraulic radius; k,is the local loss coefficient if any; y, — y, is the 
difference in hydraulic levels in a section as seen in Figure 8; and Q, the 
flow is known. By subtracting the composite roughness in the n, section. 


From Eq. 11, it follows that: 


1 
2 k 
201 y A 
n,= Y de AZ 2) (12) 
a Et) 
AZRhx /3 AZRhz /3 
Where: 
Ax = (x2— Xx1) (13) 


The roughness is composed of hydraulic concrete in the floor and 


chamfer and shotcrete in walls and vault (Figure 9). 
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chaflán 


Figure 9. Portal section with composite roughness. 


Once the roughness values n, are obtained for each section, the 
discretization of the equivalent Nikuradse roughness k, is carried out to 


obtain the real wall and floor roughness value. 


Part 2. Equivalent Nikuradse roughness estimation 
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Prandtl and Von Kármán's logarithmic law of velocities 


The logarithmic law by Prandtl and Von Kármán establish that the velocity 
profile near the wall or in the boundary layer zone can be expressed as 
follows (Schlichting, 1979): 


=In=+B (14) 


Where v represents the Reynolds-averaged longitudinal velocity; v, 
the fluid velocity in the zone of influence of the wall; y the distance to the 
wall; k=0.4 the Von Kármán constant; k, the Nikuradse equivalent 
roughness and B=8.5 correspond to the rough bottom channels. 


Therefore Equation (4) for completely rough channels is as follows: 
Y =25in2+8.5 (15) 
D, Ks 


There are three conditions for determining the constant B, which 
depends on the Reynolds number, in terms of v,k/v, where it is equal to 


5 for completely smooth walls; 5to 70 corresponds to the transition from 
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a hydraulically smooth to a completely rough regime; and greater than 
70 for completely rough flow (Schlichting, 1979) (Figure 10). 


cr ho 


------t 


A! 


02 04 : 6 5 20 22 24 26 28 30 32 
ty y 


Figure 10. Function B for roughness in terms of v,k,/v, for Nikuradse 


sand roughness (Schlichting, 1979). 


Boundary layer and displacement thickness 
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Boundary layer of a flow is the zone in which the fluid motion dominantly 
experiences the braking effect caused by the presence of the wall. The 
concept of a no-slip boundary condition is due to L. Pranatl. 

It can also be described as the necessary thickness at which the 
wall has to be displaced to make the volume lost equal to the volume lost 
due to the viscous effect of the wall (Figure 11). The vertical displacement 


is called "displacement thickness", and is defined as: 
x lo) Y) 
o =]1 (=> Jdy (16) 


Where v is the longitudinal velocity of the flow at a distance from 
the surface and the wall, measured perpendicularto the wall, and v is the 


velocity outside the boundary layer. 
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Vel. Max 


Distribución de 
velocidades en 
presencia de la pared 


h (2) 


Distribución de 
—— velocidades en 
ausencia de la pared 


Flujo de volumen perdido 


debido a la presencia de la Desplazamiento virtual de 


la pared 


ASAS 
HONDO 
ÓN 


=> DASS 


Figure 11. Interpretation of displacement thickness. 


The value of the displacement thickness 5" usually varies from one- 


tenth of the boundary layerthickness onward, according to the magnitude 
of the Reynolds number, where small viscosity ranges are large Reynolds 


numbers. 


The velocity profile in the boundary layer 
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Considering the above, the experimental methodology was proposed 
based on the Doppler velocity recordings for velocities from 10 cm up to 
2.1 m distance from the wall and with measurement segments every 20 
cm. 


First, it is necessary to verify if the recorded velocities are within 
the displacement thickness, Óa, through the following equation: 


h 
E kB Ks Int FBk=1 
ód = (1 . e a E (17) 


Where, 0d, is the displacement thickness; k = 0.4 the Von-Kármán 
constant; B = 8.5 for rough bottom channels; k, the Nikuradse equivalent 


roughness; h is the mean value of the section strain. 


Each velocity recorded by the Doppler sensor considers different 
distances along the length of the displacement thickness that generates 
different velocity profiles along with the influence of the boundary layer 
and the displacement thickness. Values that develop as logarithmic 


velocity curves are considered. 


Nikuradse equivalent roughness, ks 
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Equivalent Nikuradse roughness is obtained from Equation (5), where 


parameter Kks from logarithmic law regards: 


_L 
kk. = e (eb a) (18) 
where: 
mm _ (PY 
f =D, (E). (19) 
n=Iiní=iIn y (20) 


B defines the average velocities from the first measurement 
reading, ranging from 10 to 30 cm and afterward 20 cm every segment. 
Thus Dd, — d, refers to the velocity value obtained with the sensor at 30 cm 
(v,) and at 10 cm (0,); yn =In<3 =In y, where y is the distance in regards 


to the wall from the sensor (at 10 and 30 cm). 


And the velocity in the wall region (v,) is defined as: 


V, =04=-=—= (21) 
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Once the equivalent roughness k, is obtained, it is possible to 


calculate the value of the parameter f, known as Nikuradse's f. 


According to Aldama and Ocón (2002) in channels with uniform flow 
Sf = S, a relationship can be found between the Manning and Darcy- 
Weisbach formulas. Therefore the following expression relates the 


roughness factor (n) to the dimensionless friction factor (f): 
n= ($/89)2Rn/e (22) 


Likewise, ina flat bottom made up of sand grains, the roughness n 


is proportional to the diameter of these grains raised to 1/6 (Aldama, 8 
Ocón, 2002), which results in the following equation: 


n= kg 1?ks1/8 (23) 


Where k,, is a dimensionless constant equal to 0.129765776 
(Aldama 8: Ocón, 2012). Likewise, considering for wide channels R=h and 


clearing f we obtain: 


f =8k.,,"(ks/Ry 3 (24) 
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According to Aldama and Ocón (2002), the graphs for the friction 
factor f corresponding to circular ducts working at full or half full and wide 


channels produce values that do not exceed more than a 5 % difference, 
so the applicability ofthe above equations can be included. 


For the present study, that same procedure is being considered, 
using input values obtained from the composite roughness of Equation 
(12) and for the sections where the doppler meters are located (KO+050, 
KO+655, and K1+115). 


Likewise, those are within the hydraulic concrete roughness 
intervals between 0.012 and 0.019 and shotcrete between 0.020 and 
0.028 (Marengo 8 Arreguín, 2008). Those within the displacement 
thickness are also considered, that is, the velocity profile of each of the 
records is reviewed: only those that fit the logarithmic distribution of 
velocities along the entire length of the displacement thickness are 
chosen. 


Part 3. One-dimensional theoretical model and 
calibration of results with application to 18 empirical 
composite roughness equations 
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Once the wall roughness parameters (Equation (22)) have been obtained, 
and the composite roughness is known from the Energy equation 
(Equation (12)), then each of the 17 known empirical composite 
roughness equations (Table 3) must be applied to obtain the roughness 
parameter that best approximates to measurements. 


Table 3. Empirical equations for estimating the resistance coefficient nc 
in channels with composite roughness (Marengo 8 Arreguín, 2008). 


Ca IA O O 


Sum of component n weighted byarea Pi 
VS = ) (¿A8R5) 


ratio; orTotalshear velocityisa weighted 


sum of subarea shear velocity (1,/W) = (R¡/RO7 


Total resistance force is equal to the sum 
PyRS = XP,¡yR;S; 


of subarea resistance forces; or, n; 


weighted byvAi (MM =18,/0% 


Total discharge is the sum of s ubarea 


discharges 


Total discharge is the sum ofsubarea Q =VA = 214,320; 


2,(4,/m;) discharges (S¡/S) = (R/RJY? 


Same as HortonandEinstein's Eq. Ebut 


derivederroneously 


3| Total cross-sectional mean velocity equal 


to subarea mean velocity 
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CS O 


Total discharge es the sum of subarea Q =20Q; 


discharges (S¡/S) = (R/RI/ 


Total resistance force, F, is the sumof 


subarea resistance forces, 2Fi 


(1,JV) = (28,0% 


í 
Total shear velodtyis a weighted sum of y 9RS = y E v iS) 


(nm; P,) subarea shear velocity; or, contributing 


a (1,/V) =(R¡/R)Y? or 
componentroughnessis linearly 


proportional to the wetted perimeter 


Total resistance force, F,is the sum of 


subarea resistance forces, 2Fi 


Total resistance force equal to the sum of 
subareas resistance forces 
(V,/M = (R¡/R)Y? 


Total discharge isthe sum of subarea Q =VA= Z(V,A/) 
discharges (S¿/S) =(R/R) 


Q =VA= X(1,A4;) 
Total discharge is the sum ofsubarea 
: (S¡/S) = 1 
discharges 


R =A/P 
Q=VA=X(VA) 


(S¡/5) = 1 
Same asEq.Lwitha modified definition of 


Rde 
R 


PR3/3 AR?/3 
EPRREO  yA¡RA 
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A A 
Total shear veloaty, ./gRS ¡is a weighted gRS = MS gr) 


sum of subarea shear velocity (V,/V) =1 
f = 


Total shear velocityis a weighted sum of y9RS = y (E /ar5,) 


subarea shear velocity 
(V¿/V) = (R¡/RYP 


P 
En; PRI) Total shear velodtyis a weighted sum of y 9RS = y E y 9R5;) 


subarea shear velocity ( ol y1/2 
V¿/V) = KR¡/R 


y P a? in El Logarithmic velocity distribution over 


depth h for wide channel 225 | 1093h; 
ES in (— o) 


ki 


n =0.0342k 


This is achieved with a one-dimensional theoretical model that takes 
these equations and calibrates the logs. The one-dimensional theoretical 
model has the roughness coefficient of the shotcrete nc as a variable. In 
this case of a horizontal H2 profile, the control section is located 
downstream, and the equation that solvesitis: 


FO) = Ls Z51 +Y5 +7 A e (25) 


know (Ys) 
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a q? Q?n? Ax 
id As-1*29 A 2 


unknown  (Ys-1) 


It is a quasi-permanent model because it is analyzed with flow 
records, and the theories for a gradually varying flow are applied. 


The model starts from a known downstream value, the head is 
measured by a sensor. Once the flow is known, there is a unique 
relationship forobtaining the geometric variables such as area, perimeter, 
and hydraulic radius, that were previously defined by their corresponding 
formula. Flow, wall, and floor roughness obtained from the direct 
measurements of the Doppler meters are introduced as variables for each 
case. This scheme is applied to 17 known empirical equations to estimate 
composite roughness (Table 3). That allows for the construction of a 
hydraulic profile that will be adjusted with the profile of measured level 


values from sensors. 


To select the best adjustment equation, the concept of the Minimum 
Square Error (MSE) is used. Then composite roughness obtained in the 
measurement of Doppler sections is used through Bernoulli and Manning 
Equation (Equation (12)) and is compared with the composite roughness 
obtained with the 17 empirical equations. In both of them, wall and floor 
point roughness obtained with Doppler equipment were integrated 
(Equations (18) to (24)). With the above, the MSE obtained will allow 
definingthe best-adjusted Equation amongthe 17 known ones or defining 


a new equation. 
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Results 


Measurement of hydraulic variables 


Hydraulic levels were programmed to be measured withthe sensors every 
hour, and values from 8.72 to 10.66 m of water level were obtained. 


The velocities measured with the Doppler equipment 10 cm from 
the wall and onward, had velocity ranges from 0.35 to 2.5 m/s. However, 
there were constant communication problems with the velocity equipment 
in Tunnel 2, therefore only the information from Tunnel 1 was considered 


for the roughness study. 


The flows from Doppler "time in traverse" equipment in both tunnels 
(QT1+0QT2) ranged from 117 to 575 m/s, between October 2011 and 
February 2012. 
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Composite roughness from Bernoulli and Manning's 


equations, Equation (1) and Equation (2) 


The roughness values obtained foreach analysis section are shown in the 


following graph (Figure 12), and only those within the recommended 


ranges were considered: hydraulic concrete between 0.012 and 0.019 and 
shotcrete between 0.020 and 0.028 ( Marengo 8, Arreguín, 2008). 


Rugosidad compuesta nc 
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Figure 12. Composite roughness values were obtained from Equation 


(2) and within the recommended range (Marengo 8: Arreguín, 2008). 
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Nikuradse equivalent roughness (kx, ), Equation (8), and 
displacement thickness (d), Equation (7) 


Velocities registered from Doppler equipment in sections KO+050 and 
K1+115 were modeled. No information could be extracted from the 
equipment located at KO+655. Measurements were scheduled every hour 
from October 2011 to October 2012. Each section had 4 pieces of Doppler 
equipment (Figure 6 and Figure 7) and they were programmed with the 
other measurement equipment: the modular equipment for flow and the 


sensors for measuring hydraulic levels. 


Some complications encountered were that there was not always a 
constant measurement from all the equipment for all the time intervals. 
Also, during the rainy season there was constant damage to the 
transmission, due to the ovenloads that occurred in the instrumentation 
hut, where the SCADA systems were located. This led to continuous 
repairs throughout that time. Due to the inability to maintain the 
upkeeping, it was no longer possible to continue with the analysis of the 


data in the tunnels for more data in major conditions. 


More than 200 velocity profiles were obtained for each 


measurement interval, from 1h to 24 h, in 1 year of measurement. 


Those records belonging to a logarithmic profile were selected, 
totaling more than 300 profiles between 117 and 575 m/s of flow as 


shown in the following graphs (Figure 13, Figure 14, and Figure 15). 
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Figure 13. Velocity profiles were recorded by Doppler for different 


values of flow, from the Doppler located in Upper Left Bank Tunnel 1 
(MIST1 KO+050). 
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Figure 14. Velocity profiles were recorded by Doppler for different 
values of flow, from the Doppler located in the Upper Right Bank Tunnel 
1 (MDST1 K1+095). 
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Figure 15. Velocity profiles were recorded by doppler for different 
values of flow, from the doppler located in Lower Left Bank Tunnel 1 
(MITT1 K1+050). 


Estimates were made for different sections, and it was observed 
that the displacement thickness (0d) was about 10 % of the height water 
level value measured (Figure 16). 
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Figure 16. Review of the displacement thickness dd. 


Applying Equation (22) and the methodology described in the "Nikuradse 
equivalent roughness ks " subchapter, Nparea Wall roughness values were 
obtained (Table 4, Figure 17, Figure 18, Figure 19, and Figure 20). For 
Doppler equipment located at sections K0+050 and K1+115 of Tunnel 1 
and located at different heights: Lower Left Margin Tunnel 1 (MII T1) in h 
= 2.5 m; Tunnel 1 Upper Left Bank (MIS T1) in h = 5.5 m; Tunnel 1 
Lower Right Bank (MDI T1) in A = 4 m; Tunnel 1 Upper Right Bank (MDS 
T1)inh = 7 m. 
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Table 4. Wall roughness values (Nparea) in shotcrete for doppler at 
KO+050 and K1+115. 


Section 
Section 0+050 | Mparea | h (m) Npared 
1+115 
MII T1 |E3496 | 0.0121 z 0.0278 
E3446 | 0.0154 , 0.0174 
E3449 E 0.0139 


E3488 : 0.0134 


rugosidad, n 


80 85 90 95 100105110115120125130135140145150155160165170175180185190195200205210215220225230 
Gasto Q (m?/s) 


Figure 17. Roughness-flow graph for Doppler at 2.5 m (MIT Ti at 
KO+050). Representative roughness n = 0.016. 
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Figure 18. Roughness-flow graph for Doppler located at 5.5 m (MIS T1 
at K0+050). Representative roughness n = 0.017. 
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Figure 19. Roughness - flow graph for Doppler located at 7 m (MDS T1 
at K0+050). Representative roughness n = 0.018. 
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Figure 20. Roughness - flow graph for Doppler located at 4 m (MDI T1 
at KO+050). Representative roughness n = 0.016. 


As for water level height, it was observed that the roughness does 
not change significantly, as is the case with the flow, even when it is 
increased by 30 cm. The representative roughness is in the range of n = 


0.019 (Figure 21, Figure 22, and Figure 23). 
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Figure 21. The roughness-level graph was obtained from the doppler at 
2.5 m (MIT T1 at KO+050). Representative roughness n = 0.015. 
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Figure 22. The roughness-level graph was obtained from the doppler at 
5.5 m (MIS T1 at K0+050). Representative roughness n = 0.016. 
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Figure 23. Calculation of the equivalent roughness ks and the 


parameters f and n. 


Additionally, it was observed that the roughness varies slightly in 
the reinforced zones (K0O+160 to K0+215 and from KO+380 to K0+435). 


These zones have smooth concrete since they are structurally reinforced 
zones (Figure 24). 
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Figure 24. Location of roughness zones (shotcrete and hydraulic 
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concrete). 


Two graphs are shown: one for flows greaterthan 150 m3/s (Figure 
25), and one for flows less than 150 m3/s (Figure 26), where roughness 
changes are observed. 
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Figure 25. Composite roughness along the tunnel for a Q > 150 m/s. 


Doppler at 2.5 m height. 
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Figure 26. Composite roughness along the tunnel for a Q < 150 m/s. 


Doppler at 2.5 m height. 
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Figure 27 also shows that the flow among sections remains very 
similarbetween the first and the last sensor(section) forthe values taken 
every hour so that the calculation can be taken as a permanent flow. 
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Figure 27. Flow along the tunnel (per sensor). 


The following graphs (Figure 28, Figure 29, Figure 30, and Figure 
31) show the hydraulic profile of some of the modeled and calibrated 
outflows, the blue line represents the theoretical value, and the black dots 


the measured values. 
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Figura 29. Model calibration for Q = 217.6 m/s. 
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Figure 30. Model calibration for Q = 204.06 m/s. 
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Figure 31. Model calibration for Q = 110.97 m/s. 
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Each of the three calibrated composite roughness values was 
reviewed: nc-HG, Equation (12); Nc-measurea, Equation (22); Nc-theoretical 
(Equations Table 3), and minimum squared error was identified for all 
sections with level sensors. Those with the lowest standard deviation 
value were discretized and identified. 


In the sections along the tunnel, different roughness values were 
identified, after discretizing and calibrating them regarding the 
measurements and the calculation described before with the general 
hydraulic equation, the following values were found with the minimum 
square error (MSE), which are equivalent to an average composite 
roughness for Tunnel 1 of 0.0178 (Figure 32). 
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Figure 32. Composite roughness with MSE, per section. 
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Figure 33 shows the trend of the composite roughness values 


obtained for a range of flow rates from 130 to 210 m/s, which were 


discretized and calibrated from the theoretical model and the 
methodology proposed for the 17 equations, as well as the 
measurements. 
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Figure 33. Composite roughness (nc) conceming measured flow (Q). 


The equation with the best adjustment of the 17 known equations 


(Table 3) was identified with the letter "I", for a range of flows ranging 


from 85 to 285 m3/s: 


Equation I (15) (Chie-Yen, 2002): 
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Figure 34 shows the best adjustment composite roughness 
equations nc for each flow. 
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Figure 34. Equations of best adjustment for each calibrated flow. 


Likewise, it was possible to identify that the integration of the 
mentioned methodologies allows for a better definition of the composite 
roughness parameter in addition to the empirical equations Chie-Yen, 
2002). Then through hydraulic instrumentation, velocities can be 
measured in the areas close to the wall and, complemented with water 
height level (sensors) and flow measurements, the equivalent roughness 


parameter can be had by integrating the Logarithmic Law of velocities and 
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the Nikuradse and Prandtl-Von Kármán equations, arriving finally at the 
representative composite roughness of the section or tunnel. This 
methodology can be applied to wide channels hydraulically rough and in 
permanent regime conditions where itis possible to integrate the Darcy- 
Weisbach and Manning relationship (Aldama 8, Ocón, 2012). 


Previous experimental studies have been carried out to determine 
the composite roughness in tunnels using prototype models (Marengo € 
Arreguín, 2008). In 2007 an experimental investigation was carried out, 
where the flow was compared with four tunnel models in an arch section 
working as a channel, with different roughnesses, and under different 
runoff conditions. The prototype studied four associated roughnesses: 
acrylic, sandpaper, and plastic, which resemble hydraulic concrete, 
shotcrete, and prototype rock. The following approximations were 
obtained from the experimental results: 1) for the Acrylic-Sandpaper 
combination, the best approximation was obtained with the Ida -Engelund 
M Equation; 2) for the Acrylic-Plastic combination, the best-fit equation 
was Felkel's F; 3) finally, the author recommended that due to the 
turbulence that occurs in reality, the Felkel Equations (Equation F), Yen's 
Equation (Equation H) and Yen's Equation N should also be used (see 
Table 3). 


Conclusions 


232 
Tecnología y ciencias del agua, ISSN 2007-2422, 13(2), DOI: 10.24850/j-tyca-2022-02-04 


2022, Instituto Mexicano de Tecnología del Agua 
A Open Access bajo la licencia CC  BY-NC-SA 4.0 
Tecnología y 


Ciencias¿Agua (https://creativecommons .org/licenses/by-nc-sa/4.0/) 


The approach of the developed methodology allows for the studying of 
the composite roughness through the relationship between the contact 
surface and the fluid (boundary layer zone), applying the Prandtl-Von 
Kármán boundary layer theories using the velocity profiles obtained from 
hydraulic instrumentation using measurements very close tothe wall. This 
was appliedtothe case study of Tunnel 1 of the Grijalva River. The results 
were obtained from the Nikuradse equivalent roughness ks, and with the 
relationship between the Darcy-Weisbach and Manning equations, a wall 
roughness was obtained. Aside from this, with a general hydraulic model 
and section measurements, a representative tunnel roughness of n = 
0.0178 was obtained. It was observed that along the tunnel, the 
roughness value remains constant except for places where the surface is 
changed to smooth concrete. It was also found that Yen's Equation I 
(Chie-Yen, 2002) is the best approximation. Even though the composite 
roughness did not have a significant variationinthe range of the analyzed 
flow rates (117 to 575 m/s), it is recommended to continue with the 
measurement and to apply the methodology for higher flows since 
variations in roughness could only be observed if flow increases 
significantly. It is also recommended to continue with the water levels, 
flow, and velocity records to carry out the calibration with measured data. 


Implementing hydraulic instrumentation in real conditions in tunnels and 


channels is an important research tool for the study of composite 
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roughness in the design of these structures, which improves existing 
empirical methods. Previous studies onthese topics have been carried out 
on prototype models, however, this study additionally integrates a scale 


deviation parameter. 


The use of composite roughness, ¡.e., combined roughness of 
smooth concrete and shotcrete in tunnels, ensures more efficient 
hydraulic performance intunnels of dams, guaranteeing a better hydraulic 
discharge capacity in a contingency event during dam construction since 
it allows the integration of a longer retum period (Tr). It also maintains a 
better cost/benefit ratio by not obligating a complete recoat of all the 


sections with hydraulic concrete. 


According to Marengo (2019), in recent years great importance has 
been given to durability and failure conditions in dams and their 
temporary structures such as channels and tunnels. It is necessary to 
review these events from the technical, social and legal points of view. 
Additionally, the development of numerical methods and computational 
techniques allow for arriving at better knowledge and understanding of 
hydrological, geological, and materials aspects, among others. It is also 
necessary to implement instrumentation projects for defining the 
hydraulic and structural behavior of dams. 


The consideration of the evaluation of the failure risks togetherwith 
the composite roughness analysis methods applied in this research would 
allow for better elements for designing new projects. 
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